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Abstract 



The interaction of virtual photons is investigated using the reaction e"'"e~ — >■ e"'"e~ hadrons 
based on data taken by the OPAL experiment at e"'"e~ centre-of-mass energies ^/s^ = 189 — 
209 GeV, for > 5 GeV and at an average of 17.9 GeV^. The measured cross-sections 
are compared to predictions of the Quark Parton Model (QPM), to the Leading Order QCD 
Monte Carlo model PHOJET to the NLO prediction for the reaction e"'"e~ e"'"e~ qq, 
and to BFKL calculations. PHOJET, NLO e"'"e~ e"'"e~ qq, and QPM describe the data 
reasonably well, whereas the cross-section predicted by a Leading Order BFKL calculation 
is too large. 
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1 Introduction 



The classical way to investigate the structure of the photon at e"'"e colliders is the measure- 
ment of the process 

e+(pi)e-(p2)-e+(p;)e-(p'2)X, (1) 

proceeding via the interaction of two photons, which can be either quasi-real, 7, or virtual, 
7*. The terms in brackets represent the four-vectors of the particles as shown in Fig. |I[ 

Depending on the virtualities of the photons the scattered electrons^ may be observed in 
the detector. In the case where none of the electrons is observed (anti-tagged), the structure 
of the quasi-real photon has been studied by OPAL in terms of total cross-sections [0, jet 
production 0, and charm production 0. If only one electron is observed (single-tagged), 
the process can be described as deep-inelastic electron scattering off a quasi-real photon. 
These events have been studied by OPAL to measure the QED and QCD photon structure 
functions [§-0. If both electrons are observed (double-tagged), the dynamics of highly 
virtual photon collisions is probed. The QED structure of the interactions of two highly 
virtual photons has already been studied by OPAL In the analysis presented here, the 
investigation is extended to the measurement of the hadronic structure. The results are based 
on data recorded by the OPAL experiment at LEP in the years 1998 to 2000 at e"'"e~ centre- 
of-mass energies of ^/s^ = 189 — 209 GeV, using events where both scattered electrons are 
observed in the small- angle silicon-tungsten (SW) luminometer. The measured differential 

^Electrons and positrons are generically referred to as electrons. 
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Figure 1: The diagram corresponding to the process e"'"e e~^e hadrons. 



cross-sections are compared to the prediction of the Quark Parton Model (QPM), to a NLO 
calculation of the process e+e~ c+e~ qq, to the PHOJET Monte Carlo model [D 
and to BFKL |p calculations p!0|-p!^]. A similar analysis has been published by the L3 
Collaboration [|13[ using data taken at ^/s^ = 91 GeV and ^/s^ = 183 GeV. 



2 Theoretical framework 



In this paper double-tagged events are studied, i.e. both final state electrons are scattered 
at sufficiently large polar angles^ 6i to be observed in the detector. This corresponds to the 
situation where both radiated photons which take part in the hard scattering process, are 
highly virtual. Throughout the paper, i = 1,2 denotes quantities which are connected with 
the upper and lower vertex in Fig. |I|, respectively. 

The virtualities of the radiated photons are given by = —{pi — PiY > 0- The usual 
dimensionless variables of deep inelastic scattering are defined as: 

(li(l2 qiq2 

yi = — , y2 = — , (2) 
pm P2qi 

^The right-handed OPAL coordinate system is defined with the z axis pointing in the direction of the e~ 
beam and the x axis pointing towards the centre of the LEP ring. The polar angle 9, the azimuthal angle (p 
and the radius r are the usual spherical coordinates. 
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Qi Q\ /o\ 

2gig2 2gig2 

The e+e~ centre-of-mass energy squared is given by See = {pi + ^2)^ and the hadronic 
invariant mass squared by W"^ = (gi + ^2)^- 

The kinematical variables Q^, Ui and Xi are obtained from the four- vectors of the tagged 
electrons and the hadronic final state via: 

Ql = 2E^Ei{l- cos 9i), (4) 
= l-|icos2(^,/2), (5) 

where i?b refers to the energy of the beam electrons, and the mass rUe of the electron has 
been neglected. 

In this analysis, the hadronic invariant mass, W, is obtained from the energies. Eh, and 
momenta, ph, of final state hadrons (/i), excluding the scattered electrons: 



(7) 



For the comparison of the data to BFKL calculations the following additional kinematic 



quantity, which is a measure of the length of the gluon ladder, is defined [14 



r = ln|^^ ^In -^1 =F, (8) 



where the approximation requires W"^ ^ Qj. In the analysis presented here Y will be 
approximated by Y. 

The differential cross-section for the process of Eq. |1| in the limit ^ and for small 
values of i/i is given by |15 



d'^p'd^p' / 1 - -\ 

d^cr = ^j^^ ^ £tt [o-tt + ctlt + o^tl + ctll + -^tt cos 20 - 4rTL cos 0j 



d^p'id3p2 



£tt 0-7*7*, (9) 



E[ E'2 

/:tt, (10) 



d^p'id3p2 



E[ E'2 



where is the angle between the two scattering planes of the electrons in the photon-photon 
centre-of-mass system. The cross-sections cttt, c"tl, ctlt and ctll and the interference terms 
ttt and ttl correspond to specific helicity states of the interacting photons (T=transverse 
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and L=longitudinal). The sum of these terms is (t^*y-> cross-section for the reaction 
r^*ry* _^ hadrons. The process e"'"e~ e^e~ hadrons can be written as a product of a term 
describing the flux of the incoming photons and the cross-section for the interaction of the 
virtual photons. For any bin Ltt is derived by integrating over the phase space of the 
outgoing electrons. Here Ltt describes the flux of transversely polarised photons, and only 
depends on the four- vectors gi, g2, Pi, P2 and on the mass of the electron ||I6| . 

In the case of lepton pair production e^e~ e^e~ t^l' the cross-section is completely 
determined by QED. For hadronic final states QCD corrections have to be taken into account 
and in addition the cross-sections and interference terms cannot be completely calculated 
within the framework of perturbative QCD. For final states it is found that the 

contributions from ttt and ttl are large for Q\ ^ Q\ and small values of W'^. It cannot 
be excluded that the interference terms and contributions from longitudinal photons in the 



QCD case are also large [IC]. Therefore, only comparisons to theoretical models contain- 
ing predictions for all cross-sections and interference terms are meaningful. The cleanest 
experimental quantity which can be extracted without making further assumptions about 
the interference terms is the cross-section for the reaction e"'"e~ e^e~ hadrons as given in 
Eq. ^. Dividing by Ltt the cross-section be extracted. 

The order as QCD radiative corrections to the process e^e^ — > e+e^ qq have also been 
calculated The calculation accounts for the running of as and a and allows for the 
evaluation of differential cross-sections for several variables, but, it is valid only for massless 
quarks. The LO result is somewhat higher than the QPM prediction, which is expected due 
to the use of massless quarks. Overall the NLO corrections are predicted to be small, but 
they change e.g. the shape of the Y distribution, as can be seen from Table [^. 

Recently, much attention has been given to the BFKL pomeron ||^, especially for small 
Bjorken-x deep-inelastic electron-proton (ep) scattering at HERA. With increasing centre- 
of-mass energies squared, s, and for moderate photon virtualities, large logarithms in 1/x 
are expected to affect the rise of the cross-section. Resummation of these logarithms leads to 
the so-called BFKL evolution equation from which one derives that the total ep cross-section 
should increase as s"" with ao ~ 0.5 in Leading Order (LO), where ao denotes the pomeron 
intercept. The onset of such BFKL effects has been searched for in ep structure function and 
hadronic final state data jl^. The situation is so far still inconclusive, with the strongest 
hint for BFKL effects in the data coming from the analysis of forward vr" production. 



It has been argued [|T5|,|TS|-|2^ that e"^e~ collider offer an excellent opportunity to test the 
BFKL prediction, through a measurement of cr^*^* . For sufficiently large photon virtualities 
Qi and Ql (i.e. more than a few GeV^), this BFKL calculation can be carried out without 
non-perturbative input. If additionally ~ Ql, then the evolution in is suppressed, 
allowing for a clean test of BFKL effects. This condition is to a good approximation fulfilled 
for the data sample selected in this paper. 

The leading-log resummation for 7*7* scattering, can be represented by the so-called 
'gluon ladder diagrams', as sketched in Fig. |^. The first diagram is the 77 — >■ qq contri- 
bution and the next ones give the gluon exchange contributions. The original LO-BFKL 
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Figure 2: Diagrams contributing to the leading Inl/x approximation for the 7*7* cross- 
section. 



calculations [|T4| predicted an increase of a^*^* by a factor 20 or more compared to calcu- 
lations without BFKL effects, or with only DGLAP (Q^) evolution Since then the 
LO-BFKL calculations have been improved by including charm quark mass effects, running 
of the strong coupling constant and the contribution of the longitudinal photon polarisa- 
tion states. Recently, it has become clear that the Next-to-Leading Order (NLO) corrections 
to the BFKL equation are large and effectively reduce the value of a^. A phenomenological 
determination of the Higher Order (HO) effects was presented in Ref. |^ and the resulting 
BFKL scattering cross-sections were shown to increase by a factor 2-3 only, relative to the 
calculations without BFKL effects. Since then theoretically motivated improved higher order 
calculations have been performed [0,|T2|,p3-p6[, and give similar results. 



3 Monte Carlo generators 

Several Monte Carlo generators have been compared to the data or have been used to correct 
the data for detector effects. Most relevant to the analysis presented here are the programs 
which are used to model two-photon interactions for double-tagged, single-tagged and anti- 
tagged events. The main features of the programs are briefly described here. For further 
details the reader is referred to the original publications. An overview can be found in [|T^. 

The PHOJET event generator has been used to simulate double-tagged events. There 
are significant differences between PHOJETl.lO and the previous version PHOJET1.05 in 
several aspects of the event generation, the latter producing e.g. larger cross-sections at large 
W. In this paper we compare with the most recent and improved version of the program 
PHOJETl.lO. The program is based on the Dual Parton Model, containing both hard 
and soft processes [§. The hard processes are calculated in LO perturbative QCD, and soft 
processes are modelled based on 7p, pp and pp data assuming Regge factorisation. The 7*7* 
cross-section is obtained from the 77 cross-section by extrapolating in on the basis of the 
Generalised Vector Dominance model using the description of Ref. |2^. Events are generated 
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for both soft and hard partonic processes. A cut-off on the transverse momentum of the 
scattered partons in the photon-photon centre-of-mass system of 2.5 GeV is used to separate 
the two classes of events. For this reason the generation of events with W below 5 GeV is 
known to be incomplete. As a cross-check of the PHOJET model the PYTHIA6.130 |^ 
Monte Carlo generator, based on a recent model by Friberg and Sjostrand |^9|, is used for 
the simulation of double-tagged two-photon events. 



The general purpose Monte Carlo program HERWIG5.9-|-fct (dyn) was used to sim- 
ulate single-tagged two-photon events. This model version uses a modified transverse mo- 
mentum distribution, kt, for the quarks inside the photon, with the upper limit dynamically 
(dyn) adjusted according to the hardest scale in the event, which is of order Q^. This has 
been found ||31| to give a better description of the experimentally observed hadronic final 



states. 

The PHOJETl.lO Monte Carlo has also been used to simulate anti-tagged two-photon 
events. It is known to describe satisfactorily the OPAL anti-tagged two-photon data [|I|]. 

The QPM cross-section e"'"e~ e"'"e~qq, which corresponds to the diagram labelled Fq 
in Fig. 1^, was calculated with the GALUGA |^ program, which includes all terms from 
Eq. ^. The quark masses assumed are 0.325 GeV for uds and 1.5 GeV for c quarks. For the 
region of ly > 5 GeV considered here, the cross-section depends only weakly on the chosen 
masses, e.g. the results for u and c quarks differ only slightly. GALUGA was also used to 
calculate Ltt- 



Radiative corrections are calculated with the program BDK [Q. As for the GALUGA 
Monte Carlo, the BDK program calculates the QPM cross-section with, in addition, initial 
and final state QED radiative corrections to the scattered electrons. It has been verified 
that the non-radiative cross-sections predicted by BDK and GALUGA agree with each other. 
GALUGA has more flexibility for calculating cross-sections and is therefore used to calculate 
the Ltt factors and QPM predictions. The size of the radiative corrections depends on the 
variables used to calculate the kinematics, and also to some extent on the non-radiative 
cross-section. 

Fig. ^ compares two methods of calculating the variable Y. The first method uses the 
hadronic final state to calculate W. The second method is a hybrid that calculates W using 
yi and ?/2 obtained from the electrons, if both are positive and at least one i/i is larger than 
0.25 otherwise Y is calculated as for the first method. The non-radiative QPM cross-section 
was reweighted to the leading order PHOJET cross-section, which agrees well with the 
measured cross-section (see Section 7). Fig. |^ shows the ratio of the non-radiative (Born) to 
the full radiative, detector level measured (rad) cross-section. For the fully hadronic method 
the radiative corrections are small. However, for the electron method the corrections can be 
larger than 50% at large Y values. Obviously, measurements based on the electron kinematics 
cannot be compared with models or BFKL calculations in the region Y > 4, unless radiative 
corrections have been applied. Since the actual size of the radiative corrections also depends 
on the non-radiative cross-section itself, an iterative procedure would be required to extract 
the non-radiative cross-section. The present statistics do not permit such a procedure, so in 
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Figure 3: Radiative corrections for the process e"'"e~ — > e"'"e~ hadrons as a function of Y 
for two different methods to calculate W: a combined electron/hadronic final state method 
(electrons), and a purely hadronic final state method (hadrons), as explained in the text. 
Only statistical errors are shown. 



the analysis which follows Y has been calculated only with the hadronic variables, for which 
the corrections are much smaller. 



4 The OPAL detector 



A detailed description of the OPAL detector can be found in Ref. 34 . Here only a brief 



account of the main components relevant to the present analysis is given. 

The central tracking system is located inside a solenoid magnet which provides a uniform 
axial magnetic field of 0.435 T along the beam axis. The central tracking system consists of 
a two- layer silicon micro- vertex detector a high precision vertex drift chamber, a large 



volume jet chamber and a set of z-chambers for accurately measuring track coordinates along 
the beam direction. The transverse momenta, p^, of tracks are measured with a precision of 



ap^/pT = y'0.022 + (0.0015 ■ ptY (px in GeV). 

The central detector is surrounded in the barrel region (| cos6\ < 0.82) by a lead glass 
electromagnetic calorimeter (ECAL) and a hadronic sampling calorimeter (HCAL). Outside 
the HCAL, the detector is surrounded by muon chambers. There are similar layers of detec- 
tors in the end caps (0.81 < | cos^| < 0.98). The barrel and end cap sections of the ECAL 
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are both constructed from lead glass blocks, with a depth of 24.6 radiation lengths in the 
barrel region and more than 22 radiation lengths in the endcaps. 

The small angle region from 47 to 140 mrad around the beam pipe on both sides of 
the interaction point is covered by the forward detectors (FD) and the region from 25 to 
59 mrad by the silicon-tungsten luminometers (SW) The lower boundary of the SW 

acceptance is effectively 33 mrad due to the installation of a low-angle shield to protect the 
central detector from synchrotron radiation. 

The FD consists of cylindrical lead-scintillator calorimeters with a depth of 24 radiation 
lengths divided azimuthally into 16 segments. The electromagnetic energy resolution is 
approximately 18%/V^, where E is in GeV. 

The SW detector consists of two cylindrical small angle calorimeters encircling the beam 
pipe at approximately ±2.5 m from the interaction point. Each calorimeter is made of a stack 
of 18 tungsten plates, interleaved with 19 layers of silicon sampling wafers and mounted as 
two interlocking C-shaped modules around the LEP beam pipe. The depth of the detector 
amounts to 22 radiation lengths. Each silicon layer consists of 16 wedge-shaped silicon 
detectors. The sensitive area of the calorimeter fully covers radii between 81 and 142 mm 
from the beam axis. The energy resolution is approximately 5% on both sides. It is found 
that the energy resolution is almost constant with the energy due to energy leakage and dead 
material. In this analysis, the SW detector was used for tagging the scattered electrons from 
the process e"'"e~ — > e"'"e~ hadrons. 



5 Event selection 

The data sample used in this analysis corresponds to an integrated luminosity of 592.9 pb~^ 
accumulated by the OPAL experiment in 1998 (168.5 pb'^), 1999 (208.3 pb"^) and 2000 
(216.1 pb~^) at e+e~ centre-of-mass energies ^/s^ = 189 — 209 GeV with a luminosity 
weighted average of ^/s^ = 198 GeV. Double-tagged two-photon events were selected with 
the following set of cuts: 

1. Two electron candidates should be observed, one in each SW detector, with energies 
Ei^2 > 0.4i?b and polar angles in the range 34 < 6*1^2 < 55 mrad. The angles 61^2 are 
measured with respect to the original beam direction. On each side the calorimeter 
cluster with the highest energy is taken as the electron candidate. The energy threshold 
for the electron candidates is kept as low as possible in order to access large W values, 
where BFKL effects may become significant. 

2. In order to remove events with scattered electrons in FD or in the central electromag- 
netic calorimeter, we require that there is no single cluster in these detectors with an 
energy above 0.25ii^b- 
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At least 3 tracks {Nch > 3) have to be found in the tracking system. A track is required 
to have a minimum transverse momentum of 120 MeV, and to fulfill standard quality 
cuts as given in 



The visible invariant mass, W^is, is required to be larger than 5 GeV. It is reconstructed 
from tracks measured in the central tracking detectors and the position and energy of 
clusters measured in the electromagnetic and hadronic calorimeters, as well as in the 



forward detectors FD. A matching algorithm is used to avoid double counting of 



the particle momenta in the calorimeters and tracking chambers. 

5. To reduce the background due to beam-gas interactions, the z position of the primary 
vertex |(-2o)| is required to be less than 4 cm from the nominal interaction point. Here 
{zq) is calculated as the error weighted average of the z coordinates of all tracks at 
the point of closest approach to the origin in the r, (p plane. The standard requirement 
that the distance of any track to the origin of the z axis should be less than 30 cm was 
not applied in order to keep all possible tracks including those from off-momentum 
electrons interacting with the beam pipe. We also require that the distance of the 
primary vertex from the beam axis should be less than 0.5 cm. 

6. In order to ensure that the event is well contained in the detector and to reduce 
background from beam-gas interactions, it is required that the z component of the 
total momentum vector of the event, | Z]Pz| is less than 35 GeV, and that the total 
energy measured in the event is less than 2.2Et,. 

7. Remaining Bhabha-like events (i.e. a Bhabha event with random overlap of hadronic 
activity) are tagged using the back-to-back topology of the scattered electrons, if both 
have an energy larger than 0.7£'b. Events are rejected if the difference in radius, Ar and 
difference in azimuthal angle, A0 of the position of the two clusters in the laboratory 
reference frame are |Ar| < 0.5 cm and (vr — 0.1) < |A0| < vr rad. 



With these cuts 175 events are selected in the data. Among these we expect background 
events from two main sources. The first contribution is from e"^e~ processes containing 
electrons in the final state and the second stems from coincidences of e^e~ processes without 
electrons in the final state with off-momentum electrons from beam-gas interactions. 

The background from e"'"e~ reactions containing electrons in the final state, which amounts 
to 18.3 events, dominantly stems from the processes e+e~ —>■ e+e~rr and e~'"e~ — > e"'"e~e"'"e~. 



and was estimated using the Vermaseren Monte Carlo program ^8[. The contribution from 
other e"'"e~ background, such as e^e~ — >■ qq and other processes leading to four fermion final 
states was found to be negligible, as was the background from single-tagged two-photon 
processes with a cluster from the hadronic final state misidentified as the second electron. 

Beam-gas interactions can result in off-momentum electrons observed in the SW detec- 
tors, faking final state electrons from the process e"'"e~ e"'"e~ hadrons. This background 
was estimated using a sample of Bhabha events, selected by requiring events with two back- 
to-back electrons in the SW calorimeters, which each have an energy of more than 0.7i?b, 
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with |Ar| < 0.5 cm and (vr — 0.1) < |A0| < tt rad. Additional clusters in the SW detec- 
tors, which fulfill the criteria for electron candidates mentioned above, but do not belong 
to the Bhabha event, are counted as off-momentum electrons. The probabilities to have an 
overlapping off-momentum electron with an event coming from the interaction region are 
determined for the left {—z) / right {+z) side of the detector separately and amount to 
0.000715/0.00115 (1998), 0.00139 /0.00279 (1999) and 0.00117 /0.000671 (2000). The rela- 
tive statistical precision of these probabilities is 2-3%. The background expectation is more 
than a factor two larger in 1999 than in the other years. Using these numbers, and assuming 
that the overlap probabilities are independent between the left and right side of the detector, 
we predict that there should be 14.6 events with 'double' overlaps in the Bhabha sample, 
which agrees well with the 19 events observed. 

It is essential to check this method of estimating the background from off-momentum 
electrons on a different process. Here we used a sample of single-tagged two-photon events. 
The sample is selected with the same cuts as described for the double-tagged selection, except 
that only one scattered electron is required and cuts 6) and 7) are not applied. The event 
properties are compared with the absolute prediction for the single-tagged plus background 
events. The single-tagged events are based on the HERWIG generator with the GRV |^ 



parametrisation of the photon structure function F2, which has been shown to describe the 
single-tagged cross-section within about 10% for E > 0.7-Eb- The genuine single-tagged 
events are complemented by the physics background and by artificially created single-tagged 
events constructed by a combination of anti-tagged two-photon events generated by PHO JET 
with a cluster created by an off-momentum electron. The 6, and energy dependence of the 
electron clusters are given by the spectra of the additional clusters in Bhabha events. 

Fig. ^ shows the 1999 data, which is expected to have the largest background from off- 
momentum overlaps, compared with the prediction resulting from the sum of HERWIG and 
the background, normalised to the luminosity of the data. Results are shown for 6, and 
energy of the scattered electron normalised to the energy of the beam electrons. Fig. ^ shows 
separately for the 1998, 1999 and 2000 data samples, the missing longitudinal momentum 
and missing transverse momentum in the event, calculated including the untagged electron, 
which has been assumed to have zero transverse momentum and an energy equal to the 
energy of the beam electrons. Both distributions have been normalised to the energy of the 
beam electrons. For energies below 0.6i?b the off-momentum background clearly dominates 
and the observed angular dependences, especially in 0, clearly follow the expected shape. 

The agreement between data and prediction is very good for all variables examined, 
providing confidence that the background from overlap off-momentum beam electrons is 
under control to a level of about 10%. The off-momentum background estimate was used 
to calculate the contribution of fake double-tagged events, resulting from the overlap of one 
background cluster with a single-tagged two-photon event and the overlap of two background 
clusters with an untagged event. In total 4.3, 15.2 and 4.6 overlap events are predicted for 
the 1998, 1999 and 2000 data samples. After subtraction of all backgrounds 132.6 events 
remain, with each year contributing with 31.9±6.5 (1998), 63.5±9.3 (1999) and 37.2±7.0 
(2000) events, whereas 42.5, 53.7 and 55.5 events are predicted by PHOJET. The difference 
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between the number of events observed in the three years is of the order of twice their 
statistical error. Since the statistics of this measurement are small, the data of all years will 
be combined in the following. 

The double-tagged events are triggered by two groups of independent triggers. The first 
trigger is based on the energy deposits of the observed electrons in the SW calorimeters. 
The second group only relies on the observed tracks and clusters from the hadronic final 
state. Based on these two independent groups, the trigger efficiency of the hadronic final 
state alone, for events with two electrons with energies above 0.4ii^b, has been determined 
from the data to be (96 ±4)%. 



6 Properties of double-tagged 7^7^ events 

A PHOJET Monte Carlo sample is used to correct the data for acceptance and resolution 
effects. It is therefore essential that the shape of all important distributions is well reproduced 
by the Monte Carlo simulation. In this Section a comparison is made of data distributions 
with predictions from PHOJET. Variables calculated from the scattered electrons as well as 
variables calculated from the hadronic final state are studied. The integrated luminosity of 
the Monte Carlo sample amounts to approximately 40 times that of the data. All Monte 
Carlo distributions shown in this Section are normalised to the data luminosity. In all plots 
involving both sides of the SW detector, the sum of the distributions obtained separately for 
each side is shown. 

In Fig. ^ variables which are based on electron quantities, and the variables Xi^2, are com- 
pared with predictions of PHOJET and background estimates. All variables, the normalised 
electron energies Ei^2, the polar angles 6*1,2, the azimuthal angles 0i,2, and photon virtuali- 
ties Qi2 5 reasonably well described by the sum of the signal as predicted by PHOJET 
and the estimated background from overlaps with off-momentum electrons and other physics 
processes. Note that PHOJET does not contain any explicit effects from BFKL, which would 
show up in the region of low electron energies. Fig. ^) shows the logarithm of the ratio 
of the photon virtualities, ln{Ql/Q2), of the two photons in an event. This distribution is 
peaked around zero, indicating that the values of both photons are generally close to one 
another, which is ideal to test for BFKL effects. 

In Fig. 1^ distributions are shown which characterise the hadronic final state in double- 
tagged two-photon events: the number of tracks, A^^ch, the visible hadronic invariant mass, 
W^is, the hadronic energy, -Ehad, the variable Y as well as the sum of the longitudinal J2Pz 
and transverse J2Pt momenta. Within statistics, the agreement with PHOJET is reasonable. 
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7 Results 



The cross-section for the process e"'"e~ — > e"'"e~ hadrons has been measured in the kinematic 
region defined by the scattered electron energies Ei^2 > 0.4i?b, the polar angles in the range 
34 < < 55 mrad with respect to either beam direction, and W > 5 GeV. Differen- 
tial cross-sections are presented as a function of x, Q^, W, and the azimuthal correlation 
between the two electrons A0. Here refers to the maximum of Qf and Q2, and x is 
the corresponding value of Xi or X2- For the comparison with BFKL predictions we also 
present the differential cross-section as a function of Y. From the measurement of the dif- 
ferential cross-sections of the process e+e" — >■ e+e" hadrons we extract the cross-sections for 
the reaction 7*7* — > hadrons as a function of the variable under study, using Ltt, obtained 
separately for each bin using Monte Carlo. Technically this was done by setting cx^*^* = 1 
in GALUGA and integrating (Eq. |^) using only Ltt for each bin within the experimental 
phase space restrictions. The cross-sections for the reaction 7*7* hadrons predicted by 
the models are calculated using the same Ltt factors. For all variables except the results 
for the cross-sections of the reaction 7*7* — ;> hadrons are given at an average value of Q^, 
(Q^) = 17.9 GeV^. Due to limited statistics in the data, a simple bin-by-bin method was 
applied to correct for detector and selection inefficiencies. The efficiency. Re, and purity, -Rp, 
are defined as: 



TyDetgjHad TyDetcgiHad 

where jv^^*'^^'^'^ is the number of events which are generated in a bin and measured in the 
same bin, N^^*^ is the number of events which are generated in a bin and N^^^ is the number 
of events measured in a bin. In both definitions the terms 'generated' and 'measured' denote 
events which pass all selection cuts at the hadron or at the detector level, respectively. The 
correction factor jvHad^jyDet jg obtained by dividing purity by efficiency. For the W variable 
the purity is typically around 60% over the whole range, and the efficiency is in the range 
of 30-50%. Similar numbers are obtained for Y and x, while for the A0 and variables 
the efficiencies are around 60% and purities around 80%. The correction factor is typically 
around 1.5 and fairly constant. 

The systematic error of the measurement has been evaluated taking into account several 
contributions. All changes, except the change of the cut on Nch, are applied to the Monte 
Carlo event samples because of low statistics of the data. 



1. The error due to a possible shift of the energy scale of the SW detectors was taken into 
account by scaling the electron energy by ±1%, in accord with the uncertainty in the 
scale observed in single-tagged events, which is conservative because the differences 
observed for Bhabha events were much smaller. 

2. The uncertainty in the description of the energy scale of the central electromagnetic 
calorimeter ECAL was taken into account by varying the energy scale by ±3% [0. 
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3. To estimate the uncertainties due to the selection cuts, we have performed the following 
checks, reflecting the resolution of the variables: 

(a) The lower cut on 612 was changed from 34 mrad by ±0.4 mrad. 

(b) The lower cut on Ei^2 was changed by ±5%. 

(c) The cut on W was changed by ±0.5 GeV. 

(d) The cut on |(-2o)| was changed by ±1.0 cm. 

(e) The cut on the distance of the primary vertex from the beam axis was changed 
by ±0.1 cm. 

4. The PYTHIA Monte Carlo was used instead of PHOJET to correct the data. The 
differences between the models are generally within 10-15%, except for the high W 
region where the difference amounts to 25%. The full differences were taken as the 
errors. 

5. In modelling the region of < 5 GeV PHOJET is incomplete. This affects the 
measurement only through migrations from lower to higher values of W . To estimate 
the influence of the imperfect modelling of that region in PHOJET we have chosen the 
following method. The cross-section for events with W < b GeV have been rescaled 
such that there was a smooth transition in the W distribution around W = b GeV. 
Then the number of events generated for W <b GeV were varied by ±50%. 

6. The uncertainty in the estimation of the off- momentum background was taken to be 
10% of the background, motivated by the level of agreement observed for the single- 
tagged events. 

7. The cut on the number of charged tracks has been changed from N^h > 3 to A^^ch > 4. 
This change has been applied simultaneously to data and Monte Carlo. 

8. The 4% uncertainty in the trigger efficiency was added to the total systematic error. 

The total systematic error was obtained by adding in quadrature errors coming from the 
above checks, where the sum has been performed separately for positive and negative con- 
tributions. The main contributions to the systematic errors come from changing the Monte 
Carlo model from PHOJET to PYTHIA (—16% change in the total cross-section), changing 
the cut on A'ch (—12%), and varying the lower cuts on W (±8%) and 61^2 (±5%). The 
normalisation uncertainty due to the luminosity measurement is less than 1% and has been 
neglected. 

The total measured cross-section for the process e^e~ e^e~ hadrons in the previously 
defined phase space, is 0.35 ± 0.04 (stat) toll (sys) pb. The expected cross-section from 
PHOJET is 0.39 ± 0.02(stat) pb, while the GALUGA prediction for QPM using massive 
quarks is 0.27 ± 0.02(stat) pb, and the NLO predictions for the reaction e"'"e~ e+e~ qq 
using massless quarks is 0.35 pb. In Fig. ^we show the measured differential cross-section for 
the processes e"'"e~ e^e~ hadrons and the cross-section for 7*7* — ;> hadrons as a function 
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of X and Q^. The numerical values are given in Tables ^ and ^. PHOJET better describes 
the data at lower x values, where the QPM prediction is too low. In contrast, in the large x 
region the QPM prediction is sufficient to account for the data. 

Fig. H shows the measured cross-sections as a function of W and A0. The model pre- 
dictions indicate a slightly different shape than is observed for the data. Studies for HERA 
have shown |^0| that angular variables similar to A(j) can be sensitive to the presence of 
BFKL dynamics, but so far no calculations are available for 7*7* scattering. The data show 
that the A0 behaviour of the cross-section for the reaction e"*"e~ — > e^e~ hadrons is flat, 
while the cross-section for the process 7*7* hadrons increases from A0 = to A0 = vr. 
PHOJETl.lO does not describe the A(f) distribution, whereas QPM reproduces the shape of 
the distribution. One should remark that the earlier version, PHOJET1.05, follows the data 
in both shape and normalisation. 



In Fig. 10 we compare the measured cross-section for the processes e"''e~ — >■ e"'"e~ hadrons 
and 7*7* — > hadrons as a function of Y with the PHOJET Monte Carlo, the QPM calcula- 
tion, the NLO calculation for the reaction e"'"e~ e"'"e~ qq, and numerical BFKL calculations 



(Tables § and The BFKL predictions are shown for the LO-BFKL and a NLO 



BFKL []TT| calculation both using Y. Also shown is a (partial) HO-BFKL calculation 
using Y. All BFKL predictions are shown for Y orY > 2, except for the NLO-BFKL calcu- 
lation which has been evaluated for Y > 1. Both PHOJET and QPM using massive quarks 
describe the data equally well. Also the NLO calculation for the reaction e"'"e" e"'"e~ qq, 
evaluated for five massless quarks and using A|^^ = 0.2275 GeV, is in accord with the data. 
As can be seen from Table |I] the predicted differential cross-sections as functions of Y and 
Y are very different at small values, but get much closer at higher values. This is expected 
from the approximation made in Eq. ^. It means that at low values of Y the comparison 
of the experimental result with predictions based on Y is rather uncertain, whereas as high 
values the uncertainty from using different definitions is small. 

For all BFKL predictions shown the cross-section is significantly larger than the PHO- 
JET prediction for Y > 3, and the differences increase with increasing Y. The LO-BFKL 
calculation predicts a cross-section which is too large compared to the data. This LO-BFKL 
calculation (Bartels99) |l^ already incorporates improvements compared to the original re- 
sults |]T^ by including effects of the charm quark mass, the running of the strong coupling 
constant and contribution of longitudinal photon polarisation states. Hence BFKL effects 
as large as predicted by the LO-BFKL calculation are not in agreement with the data. BFKL 
cross-sections have been calculated to NLO (Kim99) using the BLM optimal scale 



setting. At the highest Y value the NLO-BFKL cross-section is a factor seven larger than 
the PHOJET prediction. The data lie in between these two predictions. Finally, the calcu- 
lation (Kwiecinski) contains the dominant contribution of the higher order corrections 



via the so-called consistency constraint, to all orders. Its prediction in the highest reachable 
Y range is about a factor two lower than for the NLO-BFKL calculation, and this prediction 
and PHOJET are roughly equally compatible with the data 
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8 Summary and conclusions 



A data sample collected by the OPAL experiment at LEP based on an integrated luminosity 
of 592.9 pb~^ and for e+e~ centre-of-mass energies ^/s^ = 189 — 209 GeV has been used to 
study interactions of virtual photons. Differential cross-sections for the processes e"'"e~ 
e"'"e~ hadrons are measured as functions of x, Q^, W, A0 and Y. From this, the cross-sections 
for the reaction 7*7* hadrons as a functions of these variables have been derived. Both 
PHOJETl.lO and QPM describe the data equally well for the cross-sections in x, Q^, W, 
and Y. PHOJETl.lO does not describe the A0 distribution, whereas QPM reproduces the 
shape of the distribution. Also the NLO calculation for the reaction e"'"e~ — >■ e"'"e~ qq is in 
accord with the data. 



Within current calculations BFKL effects could lead to an increase of the cross-section 
by a factor up to 20 for the largest Y values. The data rule out BFKL cross-sections which 
are as large as those predicted by LO-BFKL and NLO-BFKL calculations as presented 
in |Tl|. Calculations including dominant higher order corrections predict smaller effects in 
the LEP range and are found to be consistent with the measured cross-sections. The limited 
statistics and available Y range of the data prevent establishing or ruling out the onset of 
BFKL dynamics in this reaction. 
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range 


da /AY 


da /AY 


LO 


NLO 


LO 


NLO 


- 1 


0.071 


0.065 


0.015 


0.014 


1-2 


0.135 


0.128 


0.149 


0.138 


2-3 


0.087 


0.089 


0.111 


0.122 


3-4 


0.041 


0.047 


0.047 


0.053 


4-6 


0.010 


0.013 


0.011 


0.014 



Table 1: Predictions for the differential cross-section for the process e"'"e — >■ e"'"e qq as 
functions of Y and Y using the calculation from |]^. 
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{x) 


range 


iVev 


^ OPAL 

da; 

[pb] 


Statistical 
error 


Systematic error 


^ PHOJET 

ax 

[pb] 


up 


down 


0.06 


0.0 - 0.1 


45.6 


1.46 


0.27 


0.11 


0.26 


1.29 


0.15 


0.1 - 0.2 


50.1 


1.22 


0.20 


0.10 


0.19 


1.26 


0.26 


0.2 - 0.35 


34.9 


0.55 


0.11 


0.11 


0.20 


0.85 



iQ') 

[GeV2] 


range 
[GeV2] 




^OPAL 
[pb/GeV^] 


Statistical 
error 


Systematic error 


^ PHOJET 

[pb/GeV2] 


up 


down 


13.6 


10 - 16 


32.3 


0.016 


0.003 


0.002 


0.003 


0.018 


18.9 


16 - 22 


59.4 


0.026 


0.004 


0.003 


0.006 


0.027 


24.4 


22 - 27 


25.0 


0.013 


0.003 


0.001 


0.003 


0.020 



(W) 
[GeV] 


range 
[GeV] 


A^ev 


[pb/GeV] 


Statistical 
error 


Systematic error 


^ PHOJET 
[pb/GeV] 


up 


down 


7.2 


5 - 10 


65.1 


0.031 


0.004 


0.0053 


0.0098 


0.041 


12.4 


10 - 15 


30.6 


0.017 


0.004 


0.0012 


0.0015 


0.021 


20.6 


15 - 35 


32.5 


0.005 


0.001 


0.0004 


0.0010 


0.004 


41.5 


35 - 50 


3.4 


0.001 


0.001 


0.0001 


0.0001 


0.000 



A0 
[rad] 


range 
[rad] 


A^cv 


^OPAL 
[pb / rad] 


Statistical 
error 


Systematic error 


^ PHOJET 

[pb / rad] 


up 


down 


0.31 


0.00 - 0.63 


20.2 


0.081 


0.022 


0.010 


0.032 


0.195 


0.94 


0.63 - 1.26 


23.7 


0.098 


0.024 


0.012 


0.035 


0.163 


1.57 


1.26 - 1.89 


31.2 


0.130 


0.026 


0.013 


0.036 


0.120 


2.20 


1.89 - 2.51 


27.1 


0.117 


0.026 


0.011 


0.020 


0.083 


2.83 


2.51 - 3.14 


30.4 


0.132 


0.028 


0.009 


0.010 


0.068 



Y 


range 


A^ev 


^ OPAL 
[pb] 


Statistical 
error 


Systematic error 


^ PHOJET 
[pb] 


up 


down 


0.5 


- 1 


25.7 


0.059 


0.013 


0.029 


0.033 


0.083 


1.5 


1 - 2 


44.9 


0.105 


0.018 


0.007 


0.021 


0.141 


2.5 


2 - 3 


37.1 


0.108 


0.020 


0.009 


0.011 


0.111 


3.5 


3 - 4 


12.7 


0.040 


0.014 


0.003 


0.009 


0.040 


5.0 


4 - 6 


11.1 


0.026 


0.010 


0.003 


0.010 


0.009 



Table 2: The differential cross-section for the process e"'"e~ e+e" hadrons in the region 
£"1,2 > OAEh, 34 < 01^2 < 55 mrad and 1^ > 5 GeV, as a function of x, Q'^, W, A0 and Y. 
The average value of x, Q'^, W and the central value for A(f) and F in a bin, value of the bin 
boundaries, number of measured events in the bin after background subtraction, value of 
the differential cross-section with statistical and systematic errors as well as the differential 
cross-section predicted by PHOJETl.lO, are given. 
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{x) 


a OPAL 
[nb] 


Statistical 
error 


Systematic error 


a PHOJET 

[nb] 


a QPM 
[nb] 


10-3 


up 


down 


0.06 


4.66 


0.85 


0.33 


0.83 


4.11 


2.79 


0.0314 


0.15 


8.92 


1.42 


0.71 


1.39 


9.18 


6.23 


0.0137 


0.26 


5.76 


1.10 


1.18 


2.10 


8.95 


6.03 


0.0143 



[GeV2] 


a OPAL 
[nb] 


Statistical 
error 


Systematic error 


a PHOJET 

[nb] 


a QPM 
[nb] 


10-3 


up 


down 


13.6 


5.63 


1.23 


0.77 


1.05 


6.30 


4.61 


0.0169 


18.9 


6.37 


0.91 


0.75 


1.56 


6.66 


4.41 


0.0241 


24.4 


3.91 


0.89 


0.36 


0.98 


5.93 


4.20 


0.0165 



{W) 
[GeV] 


a OPAL 
[nb] 


Statistical 
error 


Systematic error 


a PHOJET 
[nb] 


a QPM 
[nb] 


10-3 


up 


down 


7.2 


6.79 


0.96 


1.17 


2.15 


9.03 


6.27 


0.0227 


12.4 


7.11 


1.51 


0.52 


0.63 


9.08 


5.38 


0.0117 


20.6 


5.53 


1.14 


0.38 


1.11 


4.01 


2.95 


0.0183 


41.5 


3.10 


2.06 


0.29 


0.43 


1.38 


1.27 


0.0045 



A0 
[rad] 


a OPAL 

[nb] 


Statistical 
error 


Systematic error 


a PHOJET 

[nb] 


a QPM 
[nb] 


10-3 


up 


down 


0.31 


3.18 


0.84 


0.40 


1.25 


7.64 


2.82 


0.0161 


0.94 


4.25 


1.03 


0.51 


1.53 


7.05 


3.51 


0.0145 


1.57 


6.84 


1.38 


0.68 


1.88 


6.28 


4.74 


0.0120 


2.20 


7.44 


1.67 


0.70 


1.25 


5.31 


5.94 


0.0099 


2.83 


9.56 


2.02 


0.68 


0.76 


4.90 


6.44 


0.0087 



Y 


a OPAL 
[nb] 


Statistical 
error 


Systematic error 


a PHOJET 
[nb] 


a QPM 
[nb] 


10-3 


up 


down 


0.5 


6.11 


1.34 


2.95 


3.42 


8.60 


5.37 


0.0097 


1.5 


6.84 


1.17 


0.43 


1.39 


9.14 


6.59 


0.0154 


2.5 


7.99 


1.51 


0.69 


0.83 


8.27 


4.98 


0.0135 


3.5 


3.78 


1.38 


0.32 


0.84 


3.79 


2.98 


0.0105 


5.0 


4.97 


1.82 


0.49 


1.81 


1.71 


1.40 


0.0106 



Table 3: The cross-section for the process 7*7* hadrons as a function of x, Q^, W, 
A4> and Y. The average value of x, Q^, W and the central value for A0 and F in a bin, 
value of the cross-section with statistical and systematic errors, the cross-section predicted 
by PHOJETl.lO and QPM as well as Ltt are given. 
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Figure 4: Distributions of the electron energy normalised to the energy of the beam electrons 
(a,b), electron azimuthal angle (c,d), and electron polar angle (e,f), shown separately for the 
left (L) and right (R) side of the OPAL detector, for selected single-tagged events in 1999. 
The histograms are the predictions for the single-tagged process from HERWIG, the off- 
momentum background contribution, and the background from other physics channels. 
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Figure 5: Distributions of the scaled missing longitudinal momentum (a,c,e) and the scaled 
missing transverse momentum (b,d,f) in single-tagged events for the different years. The 
meaning of the histograms are as defined in Fig. 0. 
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Figure 6: Distributions of (a) the energies of the electrons, normahsed to the energy of the 
beam electrons, (b) the virtualities of both photons, (c) the polar angles of the electrons, 
(d) the azimuthal angles of the electrons, (e) the ratio of the photon virtualities and (f) 
the X values of double-tagged two-photon events. The histograms are the predictions for 
the double-tagged two-photon process from PHOJETl.lO, the off-momentum background 
contribution, and the background from other physics channels. 
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Figure 7: Distributions of (a) the number of tracks of the hadronic final state, (b) the visible 
hadronic invariant mass, (c) the total hadronic energy (d) the variable Y, (e) the missing 
longitudinal momentum and (f ) the missing transverse momentum for double-tagged events. 
The PHOJETl.lO predictions and backgrounds contributions are as in Fig. 
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Figure 8: Cross-sections for the process e"'"e~ — >■ e^e~ hadrons in the region £'12 > 0.4i?b; 
34 < 612 < 55 mrad and W > b GeV, and for the process 7*7* hadrons, as functions 
of X for (Q^) = 17.9 GeV^ (a,b), and as functions of (c,d). Data are shown as full dots 
in the centre of the bins. The inner error bars represent the statistical errors and the outer 
error bars represent statistical and systematic errors added in quadrature. The predictions 
of PHOJETl.lO are shown as solid lines, and those of QPM as dashed lines. 
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Figure 9: Cross-sections for the process e"'"e~ — >■ e+e~hadrons in the region £"12 > OAE]^, 
34 < 61^2 < 55 mrad and W > 5 GeV, and for the process 7*7* — > hadrons for (Q^) = 17.9 
GeV^, as functions of W (a,b) and A0 (c,d). Data are shown as full dots in the centre of the 
bins. The inner error bars represent the statistical errors and the outer error bars represent 
statistical and systematic errors added in quadrature. The predictions of PHOJETl.lO are 
shown as solid lines, and those of QPM as dashed lines. 
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Figure 10: Cross-sections for the process e"'"e~ — > e^e~ hadrons in the region Ei^2 > 0.4i?b, 
34 < 6*1,2 < 55 mrad and W > 5 GeV, and the process 7*7* hadrons for (Q^) = 17.9 
GeV^, as functions of Y. Data are shown as full dots in the centre of the bins. The inner 
error bars represent the statistical errors and the outer error bars represent statistical and 
systematic errors added in quadrature. The predictions of PHOJETl.lO are shown as the 
solid lines, that of the NLO calculation of the process e"'"e~ e"'"e~ qq as dotted lines, 
and those of QPM as dashed lines. Three BFKL calculations are shown: a LO one from 
Bartels et al. (Bartels99), NLO from Kim et al. (Kim99) using Y, and the calculation from 
Kwiecinski et al., using the consistency constraint calculated for Y. 
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